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Abstract 

Canonical forms and canonical variables for inviscid flow problems are derived. In these forms 
the components of the system governed by different types of operators (elliptic and hyperbolic) 
are separated. Both the incompressible and compressible cases are analyzed and their similarities 
and differences are discussed. The canonical forms obtained are block upper triangular operator 
form in which the elliptic and non-elliptic parts reside in different blocks. The full nonlinear 
equations are treated without using any linearization process. This form enables a better analysis 
of the equations as well as better numerical treatment. These forms are the analog of the 
decomposition of the one dimensional Euler equations into characteristic directions and Riemann 
invariants. 
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1 Introduction 


In the past decade a substantial effort has been invested in numerical solutions of the Euler equations 
in one dimension. These were based on characteristic decomposition of the flow, i.e,. a decomposi- 
tion into three subsystems each of which correspond to a simple propagation along an appropriate 
line in space-time. Efficient discretization schemes were based on this representation. We refer to 
[2] for extensive reference on the subject. This approach is usually extended to multidimensional 
problems by splitting it into a sequence of one-dimensional problems. 

The above approach is limited and does not capture the real structure of the solution for 
multidimensional flow problems. Moreover, the analysis is usually based on the time-dependent 
problem although the case of interest is the steady-state problem. Consequently, an important 
feature of the system such as its mixed nature has not been given proper attention. This paper is 
concerned with multidimensional steady-state inviscid flow problems. 

An analysis which is analogous to the decomposition of the one dimensional Euler equations 
into characteristic directions and the Riemann invariants, is presented for multidimensional steady- 
state inviscid flow problems. This methodology can be applied to more general partial differential 
equations, such as the steady state and time dependent viscous flow equations. This analysis of 
the equations break up the full problem into its fundamental (irreducible) subproblems and also 
describe the interaction between them. It leads to a better understanding of the equations and to 
new type of schemes that better represent the physical behavior. 

The analysis is presented for incompressible and compressible Euler equations both both in 
two and three dimensions. It shows in a clear way the existence of quantities that propagate along 
streamlines, and other quantities whose behavior is directionally unbiased for subsonic compressible 
flows or incompressible flows. Supersonic flows, on the other hand, both in two and three dimensions 
are hyperbolic with respect to the stream direction. A change of variables is introduced in order 
to bring the system into its canonical form, and the resulting new variables are called canonical 
variables. The canonical forms are block upper triangular form, where the diagonal blocks consists 
of the basic components of the systems and off diagonal ones represent some interaction between 
the different subsystems. 

For two dimensional incompressible Euler equations the canonical variables are velocities and 
total pressure p + q 2 /2 and in three dimensions they are velocities, total pressure and normalized 
helicity, namely, (u> • V)/(V • V). Canonical variables other than the velocities propagate along 
streamlines and form the non-elliptic part of the system. The velocity components satisfy an elliptic 
system with a forcing depending on the hyperbolic parts. 

For the the two dimensional compressible Euler equations the canonical variables are velocities 
and total enthalpy, and in three dimensions a normalized helicity, namely, ( u> • pV)/(pV • pV) is 
the additional canonical variable. Quantities other than the velocities propagate along streamlines. 
The velocities are governed by an elliptic system for subsonic regimes and by a hyperbolic system in 
supersonic regimes. The structure of the canonical forms shows in a clear way the relation between 
the full problems (incompressible and compressible Euler equations) and simplified models that 
have been studied in the past, such as (div,curl) systems in two and three dimensions, and their 
compressible versions. Canonical boundary conditions are discussed in this framework. These are 
the analog of characteristic boundary conditions for one dimensional compressible flows. 
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2 On Canonical Forms for Systems of Partial Differential Equa- 
tions 

Differential equations of constant coefficient are reduced by the Fourier transform to algebraic 
equations in terms of symbols which are in general matrix polynomials in several variables. Let 
P(D) be a differential operator of the form 


P(D)= £ a -° a (2- 1 ) 

|a]<m 

where D = (Di,...,Dd), a = (aj , . . . , aj) and D a = D“ l ...D^ d is the usual multi-index notation. 
The a a are l x l real valued matrices. We want to study solutions to the problem 

P{D)u{x ) = f(x) x £ ft (2.2) 

B(D)u(x) = g(x) xedSl (2.3) 

where fi is a domain in lR d and B(D) is an appropriate boundary condition to be specified later. 

Fourier analysis in full and half space are carried out in order to analyze the interior and bound- 
ary properties of the solutions, respectively. The symbol associated with P{D) is the polynomial 

P(x,iO= £ a a (itr ( 2 . 4 ) 

|ar|<m 

where f = (£i, . . .,£<*) and £ a = . . ,£°[ J and its use transforms the constant coefficient differential 

equation into the algebraic equation 

P(iOHO = f(0- (2-5) 

For elliptic systems where the matrix polynomial satisfies 

| dct>(.«)l > C|«| 2m (2.6) 

(2.7) 

extensive research has been done with regard to smoothness properties of the solution including the 
effect of boundary conditions. The case of hyperbolic systems of differential equations, i.e._, when 
the principal part of det P(i£) has real roots for £ ^ 0, has been studied in the context of time 
dependent problems. 

For the applications discussed in this paper one is interested in analyzing steady-state problems 
which are mixed elliptic-hyperbolic. That is, part of the system is governed by hyperbolic equations 
(with respect to some direction in space) and the rest of the system by an elliptic subsystem, and 
some interaction exists between the two subsystems. 

The first step in the analysis of a mixed system is the identification of the different subsystems, 
i.e., elliptic versus hyperbolic. Next comes the identification of the variables corresponding to 
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different type of behaviors, followed by a reduction of the system to block triangular form. The 
diagonal blocks are the basic building blocks of the system. 

We discuss two types of reduction to canonical forms. In the first case the new set of variables 
are linear combination of the old ones, via constant matrices. In the other case a linear combination 
via polynomials is used. Thus, the new variables in the second method involve linear combinations 
of the old ones and their derivatives. 

C-Reducibility 

Assume that there exist invertible matrices S,T, (independent of f) such that SP(£)T is of 
block upper triangular form, i.e., 


SP(0 T = 


( Ai(0 


V 


A 2 (0 

Aa(0 


A*(0 \ 


Afc(f ) ) 


= A 


( 2 . 8 ) 


where the subsystems PihJ = 1,...,* are irreducible. 

Taking the determinant of both sides, and assuming that det(5) det(T) = 1, we obtain 


det P(0 = det Ai(0 det P 2 2 (f) . . -det P kk (0 (2.9) 

The blocks on the diagonal, correspond therefore, to the primitive building blocks of the system, 
and the off diagonal ones represent the interaction of the different basic blocks. 

Definition 1 A matrix polynomial P(0 twW be called C-reducible if there exist invertible matrices 
S,T independent of £ such that A(£) = SP(£)T is of block triangular form. 

We will refer to this form of the matrix polynomial as the canonical form. It induces a canonical 
form for the differential operator via the Fourier transform. 

Thus, we have the following theorem, 

Theorem: A necessary condition for an l x l matrix polynomial P(i£) to be C-reducible to a block 
triangular form is that its determinant admits a factorization into lower order polynomials. 

We give some examples to illustrate some important points with respect to the canonical form. 

Example I: The symbol associated with the Cauchy- Riemann equations is 


( 


*6 

*6 



( 2 . 10 ) 


Its determinant — is irreducible, therefore, this matrix is not C-reducible. 

The factorization of the determinant into polynomials of smaller degree does not guarantee that 
the corresponding matrices are C-reducible. 


Example II: The matrix 
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( 2 . 11 ) 


f *fi *6 *6 \ 

*6 *6 o I 

( *6 o y 

whose determinant 

(^K-^ + ^ + el) (2-i2) 

is a product of two lower order polynomials, is not C-reducible. This can be seen by observing that 
the eigenvector corresponding to the eigenvalue i depends on 

P-Reducibility 

For numerical applications one may allow a more general transformation in reducing a system 
to block triangular form as suggested by the following definition. 

Definition 2: A matrix polynomial P(£) is said to be P-reducible if there exist invertible matrix 
polynomials f) such that S(()P ( 0 = A(f) where A(f) is a matrix polynomial of 

block triangular form. 

The matrix given in Example is P-reducible. This can be seen from the identity 


( 1 

0 

0 \ 

( it i 


*6 \ 

( iti 

1 

0 

0 

*6 


*6 

iti 

0 = 

Kl 2 

iti 

0 

1° 

*6 

) 

v <6 

0 

it l ) 

0 

0 

iti 


Thus > using 


10 0 \ 
0 i& 

o *6 J 


s (0 = t-\o = 


1 0 0 \ 
0 it 2 *6 

o j'6 -*€* / 


we get a reduction to a block upper triangular form given by 


(2.13) 


(2.14) 


( * i 1 0 \ 

A = K| 2 tfi 0 (2.15) 

V o o if, / 

By using the operators corresponding to 5, T for the differential problem one may end up with 
other variables than the desired ones. This reduction is still important for analyzing the equation 
but the original variables may need to be solved for as well. One can define an extended system in 
which the canonical form is combined with the transformation of variable and get a block triangular 
form as 


T -i 

0 



(2.16) 
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This form which we call the extended canonical form may be required only when the matrices 
,9, T -1 depend on £. In the other case the transformation from the canonical variables to the 
original ones is done by an algebraic mapping. In the applications discussed in the next section the 
extended canonical form will be found useful for three dimensional problem. 

Assuming that S, T transform the matrix polynomial P(f) into block triangular form, we obtain 
a similar reduction for the differential equation, namely, 

(S P(D)T)(T~ l u) = Sf (2.17) 

that is, in terms of a new set of variables T _1 u we get 


A (D)v = Sf 

The boundary conditions B(D)u = g are transformed into 


(2.18) 


£ Bij{D)vj = (Qg)i (2.19) 

and we assume that for each j = 1 the operator Pjj(D) with the boundary condition Bjj(D) 

form a well posed problem. 

The simplest set of boundary conditions are those for which the boundary operator in the 
transformed variables is diagonal, and each corresponding subproblem is well posed. We refer to a 
set of such boundary conditions as canonical. 

All the arguments made above applies to linear constant coefficient problems. The variable 
coefficient case and the nonlinear case may have triangular forms for the principal part only. 

It should be mentioned that the use of which are not polynomials can reduce P(£) to 

upper triangular form. However, the elements of this form are not polynomials, therefore, do not 
correspond to differential operators. While such a representation is still useful in analysis of well 
posedness of the boundary value problem, it is not as useful in developing numerical methods for 
the problem. The representations discussed in this paper involve transformations that correspond 
to differential operators and can be used in numerical implementation of certain iterative techniques 
that are based on these forms. These issues will be discussed in a separate paper. 

3 Canonical Forms for Incompressible Flow 

In the applications we present in this and the next section, although nonlinear, a full reduction to 
triangular form is presented, without involving any linearization. 

The steady state incompressible inviscid equations in nonconservative formulation [1] are given 

by 


divV = 0 

(V • V)V + Vp = 0 


(3.1) 


which can be written in operator form as 



where Q is a diagonal matrix operator of dimension d x d whose diagonal entries are the scalar 
convection operator Q = V * V. 

Analyzing the determinant of this system reveals its structure which is necessary for analyzing 
possible boundary conditions and discretization issues. Denoting the matrix operator above by 
Li nc , and freezing the coefficients at some constant flow, Vo, one obtains 


det Lj nc (£) = — 1C| 2 ( V o • f) (3.3) 

where, £ = = £? + •• 

The determinant which is already factored into irreducible factors reveals important properties 
of the system of the incompressible Euler equations. While the operator A is elliptic, the operator 
Q is hyperbolic with respect to the stream direction. This system is therefore of a mixed type. 
Some components of the system display hyperbolic behavior, while the rest show elliptic behavior. 
The power d — 1 in the factorization of the determinant suggests that there are d — 1 quantities 
that are governed by hyperbolic subsystems. This forms a subsystem of the full equation which is 
weakly coupled to the rest. 

We proceed by identifying the hyperbolic components of this system. Using the relation 

(V . V)v = - V(V ■ V) - V x curlV (3.4) 

and defining 

P = P + ^<7 2 9 2 = V ' V (3-5) 

one obtain Crocco’s form of the equations of motion, namely, 


div(V) = 0 

— V x curlV -f VP — 0 


(3.6) 


Let (ei,e 2 ,e 3 ) be an orthonormal basis such that ej is in the direction of V, and (ei,e 2 ,e 3 ) 
form a right hand system, i.e., ei x e2 = e 3 ,e2 x e 3 = ej,e 3 x ei = e 2 

Taking the inner product of the momentum equations with V and using V • (V x curlV) = 
(V x V) * curlV — 0, we obtain 


V-VP = 0 (3.7) 

This equation implies that the total pressure is constant along streamlines. Thus, the total 
pressure P is one of the canonical variables we are looking for, and it corresponds to a factor V * V 
in the factorization of the system. For three dimensional flows there exists another quantity that 
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is governed by a a hyperbolic operator. From physical insight, a quantity that exists in three- 
dimensional and not in two dimensional flows may have to do with vorticity. The vorticity 
defined by 


a> = curlV (3.8) 

will be shown to play an important role in the decomposition of the system. 

It can be easily verified that the vorticity component in the direction of the velocity cannot be 
determined from the momentum equations. However, we can obtain the following 


- ej • ( V x u>) + ej • VP = 0, 

j — 2,3 

(3.9) 

— (ej x V) • a? + ejVP = 0, 

j = 2,3 

(3.10) 


which gives 

Using the definition of ej we get the following equations 

— q(e 2 •«) + e 3 -VP = 0 
q(e 3 • u>) + e 2 • VP = 0 


(3.11) 


which can be interpreted as equations that determine two components of the vorticity vector. 
The other component of the vorticity is obtained from 


diva? = 0 

which is a compatibility condition (div curl=0). Decomposing u> as 


(3.12) 


u> = /IV + co 1 

• ej = 0 


(3.13) 


and substituting (3.13) in (3.12) and using the continuity equation one obtain an equation for 0, 
namely, 


V-V/3 + divu>- L = 0 (3.14) 

Therefore, 0 is also governed by a hyperbolic equation. Note that while P is constant along 
streamlines, 0 admits a more complicated structure depending on the behavior of P in cross-stream 
directions. 

Summarizing, we get the following 

divV = 0 

curlV — = 0 

V- V/3 + div(co x ) = 0 

u> x = (e 2 • o> ± )e 2 -I- (e 3 • w L )e 3 (3.15) 

-q(e 2 • w 1 -) -I- e 3 • VP = 0 
g(e 3 ■ to- 1 ) + e 2 • VP = 0 
<jei • VP = 0 
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which can be written in a matrix form as 



(3.16) 


Observe that although the upper block of the system which consists of of four equations 
(div, curl) is overdetermined, its solution is guaranteed by the presence of the 0 equation, which 
states that divw = 0. Another way of looking at it is to eliminate the divw = 0 equation and to 
regard the (div, curl) system as an equation for both V and 0. In this way we obtain a 4 by 4 
system for four unknown quantities. This fixes a well known difficulty of solving (div, curl ) in three 
space dimensions. 

For two-dimensional flows 0 = 0 and the equations take a simpler form. If ei,e 2 lie in the plane 
defined by the flow then a? 1 • e 2 = 0. Thus, u can be represented by a scalar quantity which we 
denote by u>. The canonical form of the two dimensional incompressible Euler equations reduces to 


where 



(3.17) 


Do = vD x — uD y (3.18) 

These representations of the incompressible inviscid equations reveal the structure of solutions. 
That is, the total pressure P and the normalized helicity (with respect to velocity) evolve along 
streamlines, and determine the vorticity vector field. The velocity components satisfy an elliptic 
system of a (div, curl) form. 

In this representation which we call the canonical form of the incompressible Euler equations the 
elliptic part has been separated from the hyperbolic part. The variables (V,0,P) will be referred 
to as the canonical variables. 

This form of the equations suggest a natural set of boundary conditions, which we refer to as 
canonical boundary conditions for the corresponding systems. The appearance of the Laplacian 
operator A in the factorization of the determinant calls for one boundary condition at every bound- 
ary point. The term Q d ~ l implies that additional d - 1 boundary conditions are to be prescribed 
at inflow points of the boundary. From the factorization it is natural to prescribe at inflow the 
canonical variables P,0, and the extra condition to be imposed at every point is V • n where n is 
an outward normal. Note that for two-dimensional problems 0 = 0 and therefore one boundary 
condition is omitted. 
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4 Canonical Forms for Compressible Flows 

The Euler equations in conservation form [1] in terms of the variables ( V,p,H ) are 


div(pV) = 0 

div(pV <8> V + pi) = 0 (4.1) 

diw(pVH) = 0 = 0 

The analysis of the equations is simplified if one moves to the non-conservative formulation, 

( Qc Q P Qh \ ( v\ ( o\ 

pQ grad 0 P I = 0 I ( 4 - 2 ) 

0 0 pQ ) \ H ) V 0 / 

where 

Q c = pdiv + 2g^yV-(V- V) 

Qv = %Q ( 4 - 3 ) 

Qh = §fjQ 

Denoting the matrix operator above by L ey and freezing the coefficients at (Vo,po?#o)> the 
determinant of the symbol gives 


det L e (0 = Po +2 (V o • 0"((Vo • O 2 ~ cgiei 2 ) ( 4-4 ) 

where Co is the speed of sound and can be evaluated using the relations (where subscripts are 
omitted), 


5 = p/p 7 


H = 1 + I ? 2 

= tf - w 

c 2 = 7 p/p 


(4.5) 


The factor V 0 • £, appearing n the determinant here is the same one as for the incompressible 
equation, suggesting that there are d quantities that propagate along streamlines. The rest of the 
system is governed by a system equivalent to the well known full potential operator Q 2 — c 2 A whose 
character depends on the Mach number M, 


M 2 = q 2 /c 2 (4.6) 

leading to an elliptic equation for M < 1 and hyperbolic (in the stream direction) for M > 1. As 
before we use Crocco’s form of the equation for analysis purposes. 


Q c V = 0 

— pV x — pTVs + pV H = 0 
div(pV/7) = 0. 


(4.7) 
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By subtracting H times the continuity equation from the energy equation one obtains 


pV . Vff = 0 


(4.8) 


Taking the scalar product of the momentum equation with V leads to 


pV -VH - pTV -Vs = 0 (4.9) 

Thus, two of the quantities that propagate along streamlines are the total enthalpy ff and the 
entropy s. As the determinant suggest there is one more quantity in three dimensional flow that 
propagate along streamlines. 

To find that quantity we follow a similar path as in the incompressible case by introducing the 
vorticity and its decomposition as 


u> = curlV 

w = /3pV + w x (4-10) 

w x • V = 0 

The equation for /3 is obtained the same way as in the incompressible case, i.e., using diva? = 0, 
giving 


pV • V/? + divw x = 0 


(4.11) 


The quantity o> x satisfies 


-pq(e 2 oj x ) - pTe 3 • Vs + pe 3 • Vff = 0 
pq(e 3 ■ a3 x ) - pTe 2 ■ Vs + pe 2 • Vff = 0 


in complete analogy with the incompressible case. 

The continuity equation is simplified by subtracting from it appropriate multiples of the energy 
and the entropy equations, given Q C V = 0. 

Summarizing we get 


Q c V = 0 

curlF — u> = 0 

pV • V/3 + div(u> x ) = 0 

u> x - (e 2 • o» x )e 2 - (e 3 • u> x )e 3 = 0 

-pq(e 2 -w 1 ) - pTe 3 - Vs + pe 3 - Vff = 0 

pq(e 3 • u> x ) - pTe 2 ■ Vs + pe 2 • Vff = 0 

-pTV • Vs = 0 

pV ■ Vff = 0 


or in matrix form 
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I Q 


—curl 


V 


pqel 


e2 e3 


pV • V dive2 dive3 
0 -q 0 

0 0 q 


0 

-pTe 3 • V 
-pTey • V 


-pTV ■ V 


0 

pe 3 • V 
pe 2 -V 


0 


VV / 



\ 


V 


0 


U-2 


w 3 


S 


* / 


= 0 


(4.14) 


The nature of the subsystem 


Q c V = 0 (4.15) 

curlV = w (4-16) 

can be studied by taking the gradient of the first equation and subtracting from it the curl of the 
second equation. Using the relation 


V 2 V = graddivV — curlcurlV (4.17) 

and 

V.gradp = j^V.(V.VV) (4.18) 

we get 

V 2 V + grad(i^^V • (V • VV)) = - curlu, (4.19) 

This equation is elliptic for subsonic flow and hyperbolic with respect to the stream direction for 
supersonic flows. 

In two dimensions (5 = 0 and the canonical form of the Euler equation reduces to 


where 


D\ D 2 

-Dy D x 

1 o 
o o 
o o 


( U > 
V 

0 0 

q 7 ( 7 - 1 ) D * g D 3 


u 

0 0 

0 TpQ 0 


s 

0 0 

0 0 pQ ) 


V H / 


/ 0 \ 
0 
0 

VO/ 


Di = p/c 2 ((c 2 - u 2 )D x - uvDy) 
D -2 = p/c 2 ((c 2 - v 2 )D y - uvD x ) 
D 3 = vD x - uD y 


(4.20) 


(4.21) 
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From the canonical form one can easily see that a natural set of boundary conditions for the 
steady state Euler equations is the following. 

At subsonic points one condition have to be specified on account of the elliptic part. This is 
the compressible analog of the (div,curl) system for the incompressible case. This condition can be 
pV • n where n is the outward normal the boundary. A compatibility condition has to hold, that ^ 

is, the integral of that quantity around the boundary has to vanish, similar to the incompressible 
case. At inflow boundary points d more conditions has to be specified and the natural choice is 
(#,$,/?), as these quantities are carried into the domain by the convection operator Q . At outflow 
boundary points no condition is required for the operators Q. 

At supersonic points the system is purely hyperbolic and require at the inflow d + 2 conditions, 
and no conditions at outflow. In that case it is natural to specify the quantities (V, s), rather 

than using (3 for one of them. 

5 Conclusion 

Canonical forms for inviscid flow problems have been derived. In these forms the different types 
of subsystems have been separated using a set of new variables, namely, the canonical variables. 

These forms are upper triangular operator forms in which the building blocks of the systems reside 
on the diagonal and the interaction of the different type of subsystems is represented by the off 
diagonal blocks. This decomposition of systems of partial differential equations allows a better 
insight into the structure of solutions which is essential in constructing numerical solutions. 

The forms described here suggest new iterative solutions of the inviscid equations, both in two 
and three dimensions. These will be discussed elsewhere. 

The different approximations used in fluid dynamics over the years are clearly seen in the canon- 
ical form. Assuming the flow is of constant total pressure P one obtain a (div, cur/) system, which 
in two dimensions is nothing but the well known Cauchy-Riemann equations. The three dimen- 
sional implementation of a (div, curl) system was not popular since the overdetermined system 
needs to be discretized carefully. As the canonical form suggests, the introduction of /? into that ! 

overdetermined system fixes the major numerical difficulties related to existence of solutions. 

For the compressible case, by assuming that the total enthalpy H and entropy s are. constant, 
one obtains a system which is equivalent to the full potential equation. This is the upper left block 
in the canonical form. ' 
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